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Abstract 

We report on and analyze the textural transformations in islands, thicker circular domains, 
floating in very thin free standing chiral Smectic C* liquid crystal films. As an island is growing, 
an initial pure bend texture of the c-director changes into a reversing spiral at a critical size. 
Another distinct spiral texture is induced by changing the boundary condition at the central point 
defect in the island. To understand these transformations from a pure bend island, a linear stability 
analysis of the c-director free energy is developed, which predicts a state diagram for the island. 
Our observations are consistent with the theoretical phase diagram. 
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I. INTRODUCTION 



Defects and singularities of the director field in two dimensional (2D) ordered molecular 
systems with in-plane orientational order cause fascinating textures, easily visible by means 
of polarized-light microscopy. In this paper we report a detailed experimental and theoretical 
study of a set of simple textural transformations in a system formed by blowing smoke over 
a thin free standing film of chiral smectic C (SmC*) liquid crystal. Sub-micron smoke 
particles nucleate islands, circular regions of added smectic layers, which grow to a certain 
equilibrium size, with the smoke particle as a point disclination at the center of the island. 
The island maintains azimuthal symmetry during growth, initially with the c-director field, 
a unit vector c, defined by the projection of the tilted long molecular axis onto the layer, 
tangentially oriented. We find that this initially stable pure bend texture can transform to 
two other equilibrium spiral splay-bend textures, either as the islands grow, or when they 
are momentarily perturbed by external forces. To understand these transformations among 
stable or metastable textures, we develop a linear stability analysis for this general class of 
textures, which predicts a state diagram, as a function of island radius, that agrees well with 
our observations. 

Important precursors to this work include an initial report of these phenomena^], the 
report of a reversing spiral texture first discovered in domains at the air/water interface when 



the SmC* film is both polar and chiral and study of the boojum, a point disclination 
occurring at the edge of a circular domain 0,0]. The theoretical analysis of the stability of 
textures around topological defects has been investigated for Langmuir monolayers with a 
tilted phase and for free standing Sm C films 0, 0, 0,0] ■ Especially important, K. K. Loh et 
a/.yi pointed out that there can be more than one stable or metastable texture for a point 
defect centered in a circular domain, and described the reversing spiral texture of a circular 
domain in a chiral system, which can be energetically stable or metastable. The remarkable 



growth of islands to an equilibrium size, once they are nuc 
chirality and their resulting polar (ferroelectric) symmetry. 



eated, is a consequence of their 
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II. TEXTURAL TRANSFORMATIONS OF ISLANDS 



Samples are prepared by drawing a small amount of smectic C* liquid crystal material 
across a 6 mm diameter hole in a thin metal sheet. Materials were prepared from the CS 
series of mixtures available from Chisso Co. Typical films as drawn were on the order of 10 
smectic layers thick. [10] To create islands, smoke from burning paper was wafted across the 
film. Watching the sample using a polarizing microscope with crossed polarizers, one could 
observe the growth of islands nucleated by some of the smoke particles. The islands were 
brighter than the very dark background film, their brightness being proportional to their 
thickness. Once nucleated, islands grew in area but not in thickness. As the islands grew, 
one could observe their internal texture due to their birefringence. Although sometimes the 
smoke particle was the center of a point defect at the periphery of the island, a so-called 
boojum, most often the sub- visible smoke particle was centered in the island as a -|-1 point 
disclination, as shown in Fig. ^ 

There are four possible textures for the intensity pattern shown in Fig. ^ For a SmC* film 
of thickness d with angle if) between c and either of the polarizers, with crossed polarizers, 
the intensity of transmitted light is 



where A is a wavelength of the incident light and An is a birefringence for light propagating 
normal to the film. 

We neglect the effects of the twist of c through the film thickness, since the islands are 
thin compared to the helical pitch of the spontaneous twist produced by molecular chirality. 
The helical pitch in the materials studied was in the range of 3 to 15 micrometers, while 
the islands were in the range of thickness of 0.06 to 0.15 micrometers. |10] If the c director is 
parallel to either of the polarizers, ■?/' = 0, vr or ±7r/2, then / = 0. At ip = ±7r/4, ±37r/4, / 
is a maximum in Eq. (^. Hence, the four orientations of c indicated in Fig. ^ are possible. 

To distinguish among these cases, the lower edge of the sample holder is lifted producing 
a small tilt with respect to the microscope stage. After the floating islands move to the edge 
of the film due to gravity, their appearance changes to that shown in Fig. |21 The increased 
intensity on the right half of the islands indicates that their c-director texture is tangential 
counterclockwise, since in that case, the added tilt increases the molecular tilt, relative to 
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FIG. 1: Tangential or radial? An island observed with crossed polarizers and four possible c 
director configurations, tangential (clockwise or counterclockwise) or radial (outward or inward). 



FIG. 2: Islands observed when the lower edge of the sample holder is raised with respect to the 
microscope stage by about 10°. The texture is tangential counterclockwise. 

the propagation direction of the light, and thus the apparent birefringence, on the right, and 
decreases it on the left of the center. All the +1 central defect islands, as shown in Fig. d 
have this "pure bend" texture throughout the island. 

To explore the stability of pure bend islands, which are apparently at equilibrium, we 
blew on the film with a small jet of gas. The resulting swirling often induced two kinds of 
distinct changes from a pure bend texture. In all cases the boundary condition at the outer 
edge of the island remained tangential counterclockwise. However, the boundary condition 
at the smoke particle at the core of the point defect would (a) sometimes change from 
tangential to approximately radial, or (b) remain tangential. Right after blowing, there 
was often a very tightly wound spiral in the texture, which slowly relaxed to equilibrium, as 
shown in Fig. |21 Over the course of minutes, as the c-director relaxes, in case (a) this texture 
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FIG. 3: The textures after blowing on islands with pure bend texture, (a) A pure bend island 
switches to a 'simple spiral' with radial c-director at the core, (b) A pure bend island exhibits a 
transient (unstable) reversing spiral texture, when the core boundary condition is unchanged. 

equilibrates to a final state with roughly ±7r/2 rotation from radial at the core to tangential 
at the outer boundary (Bottom of Fig. IHl^a)). We call this spiral texture a "simple spiral". 
We have observed the simple spirals with radial inward or outward boundary conditions at 
the core, and right or left handed spirals. Neither the sense of the spiral nor the sign of 
the radial boundary condition is meaningful, since observing the sample from the other side 
would clearly reverse both these characteristics. There is a size-dependence in this textural 
transformation; pure bend islands of a small size are hardly ever transformed into simple 
spirals. This indicates that there is only a small interaction energy between the molecular 
orientation and the smoke particle, since the pure bend texture is being stabilized by the 
torque transmitted from the outer island boundary to the defect core. 

However the other spiral texture resulting from blowing, case (b), which maintains its 
tangential boundary condition at the core, which we call a 'reversing spiral', goes back to 
the pure bend texture, as shown in Fig.lHfb), when it started from an apparently stable pure 
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FIG. 4: Transformation of a pure bend island into a reversing spiral island. It took around 30 min. 
bend texture. 

Remarkably, as shown in Fig. |3J as a stable pure bend island grows, at some critical size 
the c-director texture spontaneously starts to evolve into a reversing spiral. This transition 
is driven by a competition between the bend and splay elastic energies. If the splay elastic 
constant is less than the bend constant, this enables the transformation from a pure bend 
texture to a reversing spiral in a large enough island. The time it takes to reach equilibrium 
depends on the thickness of the island. The thinner the island, the less time. It took 30 
minutes for the pure bend island in Fig. |3] to reach equilibrium. 



III. LINEAR STABILITY ANALYSIS 



To explore the stability of the texture in the islands, we employ linear stability analysis. 
To find the equilibrium configuration of c in the islands, we should minimize the following 
free energy consisting of an integral of the elastic energy density over the area of the island, 
and the edge energy at the sample boundaries. 

F = j [KsiV ■ cf + Kb(V X cY]dA + j a{(P)dl. (2) 

Ks and Kb are the two dimensional splay and bend curvature elastic constants, respectively. 
is defined as the angle between c and the outward radial vector f. is the anisotropic 
energy per unit length of edge. Note that it is a periodic function, cr(0 — 27r) = (t(0). 
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Generally it is written as 

cr(0) = (To + y^[an cos(n0) + c„sin(n0)], (3) 



oo 



n=l 



in which ctq is an isotropic term In the areal free energy density, the possible 

linear terms, V • c and V x c, are converted to line energies by the divergence theorem 
and Stoke's theorem. Therefore the first order terms aicos0 and Cisin0 in Eq. (jH)) arise 
from a spontaneous splay(V • c) and a spontaneous bend(V x c) respectively. However in 
a free-standing film, which is physically the same when viewed from the top or bottom, no 
spontaneous splay of c is possible; for films floating on a liquid this term is allowed j2| by the 
polar symmetry of the sample. So ai must vanish, ci arises directly from the spontaneous 
bending (in our case, left turning) produced by the molecular chirality in the SmC*, and it is 
directly responsible for the minimum energy textures having a counterclockwise tangential 
boundary condition at the outer edge of all the islands we observed. We include in cr(0) 
only the lowest terms for the line energy, 

cr(0) = ctq + ^2 cos 20 + ci sin 0. (4) 

For the islands we have created by a heterogeneous nucleation with a strong anchoring 
tangential boundary condition at r = i? (outer radius of island) and a weak anchoring 
boundary at r = e (radius of the particle at the center), the line energy is 

y a{cf>)dl = 2nRa{'^) + 27Tea{(f){e)). (5) 

The first term on the right side is the line energy for the outer boundary and the second is 
for the inner boundary. The total free energy for a island in Eq. Q can be rewritten as the 
sum of the areal elastic free energy and the line energy of Eq. (0) in terms of x = ln(r/e). 

F = 7cK (1 + /i cos 20)1--) +2/isin20-- 

+ (1 -/i cos 20)1 d2; + 27re(T(0(O)), (6) 



where, 



Xq = ln{R/e), 

^_Ks + K, 

K = , 

2 

_ -K, + K, 



The constant term 2tiR(j{ti /2) in the total free energy is ignored, because it has no effect 
on Here we assume that the radial configuration of the c director, = or = vr, in 

the areal free energy is a state for minimum energy, based on our observation that c tends 
to rotate toward a radial orientation in large enough islands. Hence Kg < Kf,, which means 
fi ranges from to less than 1. 

The equilibrium condition for from the Euler-Lagrange equation derived from Eq. (jH)) 

is 



[I + fi cos 20) 



/i sm z(p — 

-dx 



The boundary condition at x = obtained from (^^^ 
K 



■ /i sm 

9f 



-- 0. 

is 



(7) 



[I + fj, cos 20) ( ^ ) + /i sin 20 
^ dx 



x=0 

da 



x=0 



x=0 



(8) 



where / and t] are the areal and edge free energy densities. 

We found 0(x) numerically, using strong anchoring at the outer boundary, 0(xo) = 7r/2, 
as required by our observation. However, depending on the choice of parameters, there might 
be one or more numerical solutions. Linear stability analysis is used to determine if these 
solutions define local minima for the free energy. For the analysis, we need to manipulate 
the total energy in Eq. (jHl). Defining 



A(0) = 1 +/icos20, 
5(0) = 2/isin20, 

C(0) = 1 — |UCOS20, 

Eq. (jni) can be rewritten as 

F = nK Pt/xU(0)f^)' + 5(0)^ + C(0) 
Jq L \dx/ dx 

+ 27rea(0(O)). 
Let's start with the second term in the integral. 

TiK / B{(p)-j-dx = ttK B{(p)d(t) 
Jo dx ^0(0) 

= 7rfryu(cos 20(0) — cos 20(xo)) 
= 7ri^:/i(cos20(O) + 1). 



(9) 
(10) 

(11) 



(12) 



(13) 
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We rewrite Eq. (fT^ . 



ttK 



XO 



dx 



(14) 



where |7ri^'S'((^(0)) = 27recr((/)(0)) + Trii'/i cos(2(/)(0)) + nKfi. We next incorporate the hne 
energy into the integral and symmetrize the problem around x = 0, obtaining 



XO 



M<l>){^) +C(0) + 5(a;)5(0(a;)) 



dx. 



(15) 



in which 6{x) is a Dirac delta function. 

Now we expand the above free energy to second order in a small fluctuation (f = (p — (po, 
in which (pQ is a solution that makes the total free energy F an extremum. The boundary 
conditions are ip{—xo) = (p{xo) = due to the strong anchoring. The expanded free energy 
6F is 

ttK 



6F 



^0 f r 1 



X 



+ 



dx 



I dip 



dx dx 



dx) 



C{ct>o) + C\ct>o)^ + -C"{<P,)^^ 



(16) 



+5{x) ( S{ct>o) + S'{cPo)v + i^S"{ct>o)v'' 



dx. 



where the notation ' and " denote first and second derivatives with respect to (p. 

Since we are expanding around an extremum solution, the first functional derivative 
5Fi = 0. As a result, the lowest term is the second order derivative of the free energy, 

TiK r° rl 



SFo 



-A"( 



90j 



+6{x)-S"i<Po)ip' 



dx 
1 

dx ^ dx 2 
dx. 



di^Y 
dx / 



+2A'(0o)#^^^ + ;C"(0o)^^ 



(17) 



We next convert the second derivative terms to a symmetric form by integration by parts. 



2A'( 



d(po dip 

-d^^'Tx 



d(podp> , nC^V 

-A (do)- — —(p - A{(po)-r^(p, 



dx dx 



dx"^ 



-A"i<Po) 



^o^^ 2 



dx 



-A'( 



d\ 



dx"^ 



2 
ip . 



(18) 



After substituting Eq. (fTHj) into Eq. (fTTj). we obtain 

TiK n r / , \ 



2 \ dx J dx'^ 

+ I C^"(0o) + i 5(x)S"(0o(O))l ifdx. (19) 



(5Fo has the form of 



irK r° 

SF2 = — [ipL^]dx, (20) 

where L is the operator ~"^(^(0o)^) + V{x) and l^(x) is completely determined by the 
solution (poix)- This symmetrized second derivative term in L is now Hermitian or self ad- 
joint, as is V{x). The operator L has an infinite set of eigenfunctions which form a complete 
orthogonal set of functions for constructing all suitable functions satisfying the boundary 
conditions on the integral. 12] The eigenfunction with the lowest eigenvalue does not cross 
zero between the ±Xo end points, while sequentially, each higher eigenvalue corresponds 
to an eigenfunction with one more zero crossing. Let's expand a test function ip in the 
eigenfunctions of L, = ^aj(y9j, in which is normalized so that 



^tdx = 1, (21) 



XO 

and each eigenfunction ipi has an eigenvalue Cj. Then 

= — y ^L^dx = ^ E aid (22) 

If the lowest eigenvalue Ci > 0, 5F2 is positive, and </)o minimizes the free energy with 
respect to small variations. If Ci < 0, then ipi clearly lowers the free energy and 0o is not 
a minimizing function. To test a solution 0o we find the eigenfunction of the operator L, 
which does not cross zero between the endpoints, and determine its eigenvalue. A negative 
eigenvalue indicates instability. 



IV. RESULTS 



We have described the procedures to obtain the equilibrium textures and to examine 
their stability against infinitesimal fluctuations to understand what we observe, which is 
that there can be more than one stable or metastable solution, for the texture of an island. 
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FIG. 5: A plot showing the minima of the free energy of an island as a function of the initial 
slope of (j) at the outer boundary for the following parameters, K = 1, fi = 0.13, eci = 1.0, and 
ea2 = 0.6. At initial slope 0, the free energy is for a pure bend island. 

To find minima of the free energy, we numerically solve the Euler-Lagrange equation Eq. (0) 
with the boundary conditions, Eq. (jSJ at the core and strong anchoring = 7r/2 at the outer 
boundary. We adopt the shooting method, in which starting from the outer boundary with 
= 7r/2, we vary the slope of 0, and integrate to the inner boundary. The discrepancy from 
the desired inner boundary condition is used to adjust the initial slope, repeating until the 
solution is obtained. 

In Fig. we show the free energy of functions as the slope at the outer boundary is 
varied. The local minima represent solutions of the Euler-Lagrange equation that satisfy 
the boundary condition at the inner boundary. For the four island sizes shown, compare 
the free energy of the pure bend island, with initial slope zero, to that of the other local 
minima. The global minimum at large initial slope represents a simple spiral, while the 
local minimum that appears for island size just above ln(i?/e) = 5.5 is a reversing spiral. 
Actually many local minima can exist at larger initial slope, but they have much higher free 
energy than the free energy of the pure bend island. The linear stability analysis of the pure 
bend texture shows that once the reversing spiral solution appears, the pure bend texture 
is unstable. 

An example of a stable simple spiral is illustrated in Fig. IHl In Fig. Efa), varies 
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FIG. 6: Stable solution for a simple spiral. K=l, ;U=0.13, eci=1.0, and ea2=0.6 at ln(i?/e)=6.0. 
(b) Texture of the solution as seen between crossed polarizers, (c) A simple spiral we observed, 
after blowing on the sample. 

monotonically from the core to the outer boundary. This is an illustration of the lowest 
energy configuration for /i=0.13, K=l, eci=1.0, and ea2=0.6 at ln(i?/e)=6.0. Fig. IHl^b) 
represents the appearance of the island, as seen between crossed polarizers, for which (j){x) is 
shown in Fig. IHl^a). The island in Fig. IHfc) is a typical simple spiral we created by blowing. 
Note that in the case calculated here, the stable boundary condition at the core is not radial 
but almost reversed tangential, but the part of the solution with values of < is so close 
to the core that it is essentially invisible at optical resolution, since the radius of the island 
is only on the order of 20 micrometers. The important qualitative observation is that 
increases monotonically from the core. In fact we do not know the precise core size, just 
that it is sub- visible. 

In the reversing spiral seen in Fig. (p rotates in one direction from the inner boundary 
until reaching a minimum angle, and then rotates in the other direction to match the outer 
boundary condition. This configuration does not represent an absolute minimum of the free 
energy, but is locally stable for /i=0.6, K=l, ln(i?/e)=4.5, eci=1.0, and ea2=0.5. Fig. [7|^b) 
shows the appearance of the island as seen between crossed polarizers, for the solution (j){x) 
shown in Fig. [7[a). Fig. [7|^c) is the texture of an observed reversing spiral. Again, the 
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FIG. 7: Stable solution for a reversing spiral. K=l, /x=0.6, eci=1.0, and ea2=0.5 at ln(i?/e)=4.5. 
(b) The texture seen between crossed polarizers, (c) A typical reversing spiral island. 

detailed texture near the core is unresolved, but qualitatively it is clear that the spiral 
reverses curvature at very small radius. 

Fig. IHl is a state diagram for an island, as a function of its radius, obtained by numerical 
calculation, where the simulation parameters are K=l, /i=0.13, eci = 1.0, and ea2=0.6. At 
small radius, the energy for a simple spiral is much higher than the energy of a pure bend 
texture. This explains why we can't observe the simple spirals in islands of small size. 
Above a certain size, the pure bend island is only metastable, but there is a sizable energy 
barrier for transition to the simple spiral. However, above a second critical size, the pure 
bend island becomes unstable relative to a reversing spiral with a lower energy, without any 
energy barrier. The transition size is sensitive to Kg and K^. If Kg and Ki, are the same, 
the size of the transition diverges. The bigger fi, the smaller the size for this second order 
transition. 

For island radius larger than the radius for the crossing of the free energies of the pure 
bend texture and the simple spiral, around ln(i?/e) = 2 in this case, the pure bend texture 
can be transformed into a simple spiral as a first order transition. We think blowing on the 
sample makes this possible in our experiments, because the highly curled transient spiral 
state it produces provides the torque needed to change the inner boundary condition. The 
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FIG. 8: State diagram for an island. K=l, /i=0.13, eci=1.0, and ea2=0.6. 



spiral then relaxes into the stable simple spiral state. 



V. DISCUSSION 



We have observed the textures of circular islands of a ferroelectric smectic C* liquid 
crystal, especially studying first and second order transformations of textures as a function 
of island size. First our basic observations confirm the fact that the chirality of the smectic 
C* phase induces a preferred bend curvature in the plane of the smectic layers, leading in the 
case of the materials we studies, to a general left turning bend texture as their ground state. 
Second, using an elastic model in which the difference of elastic constants, Ki, — Ks, provides 
a driving force, and the anchoring of the c-director at the core of the central disclination is 
weak, we construct a plausible state diagram that agrees with our observations. The circular 
geometry of the island, combined with the strong tangential anchoring of the c-director at 
its outer boundary, demand a curved texture, which makes the difference between the splay 
and bend curvature elastic constants crucial. The continuous transition from the pure bend 
island to the reversing spiral is completely analogous mathematically to the behavior of 
a plane parallel sheet of nematic with tangential boundary conditions, in a perpendicular 
magnetic field. As the sheet grows thicker, at a critical thickness, a continuous Frederiks 
transition occurs, in which the director rotates toward the external field direction. Here, 
that effective field direction is radial, converting bend into splay, and lowering the energy. 
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In studying different materials, using the textural transformation from pure bend to 
reversing spiral as a diagnostic, we saw a strong correlation between the spontaneous elec- 
trical polarization of the ferroelectric phase and the apparent difference between the bend 
and splay elastic constants; high polarization materials matched our model of a material 
with a high bend elastic constant. This result correlates well with the fact that bend cur- 
vature of the c-director produces divergence of the spontaneous polarization, resulting in 
space charges which interact to increase the free energy of the bent state. In low frequency, 
or quasi-static, conditions, the resulting space charge is partially screened by free charges in 
the material, reducing the otherwise long range electrostatic interactions to an effective local 
contribution to the free energy, which can be expressed as an increase in the effective bend 
elastic constant. [13] The fact that our state diagram and the textures we calculate from the 
elastic model agree well with our observations leads us to hypothesize that we have success- 
fully approximated the electrostatic interactions and energy by using a large bend elastic 
constant Ki, to account for the behavior of high polarization materials. The high value of K^, 
pushes the critical island size for the appearance of the reversing spiral to very small values, 
as one would expect qualitatively, and as we do observe for high polarization materials. To 
explore this hypothesis further, we made independent studies of the elastic constants and 
the interactions of bend distortions with free charges in free standing SmC* films, com- 
paring high frequency and quasi-static behavior. To accomplish this, we performed light 
scattering experiments on c-director fluctuation dynamics, with electric field quenching of 
fluctuations, with analysis including both bend-induced polarization charge and conduction 
charge screening. They will be reported elsewhere. 
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